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ABSTRACT 

Wc have analysed 10 months of public data from the Kepler space mission to measure 
rotation periods of main-sequence stars with masses between 0.3 and 0.55 Mq. To 
derive the rotational period we introduce the autocorrelation function and show that 
it is robust against phase and amplitude modulation and residual instrumental sys- 
tematics. Of the 2483 stars examined, we detected rotation periods in 1570 (63.2%), 
representing an increase of a factor ~30 in the number of rotation period determi- 
nation for field M-dwarfs. The periods range from 0.37-69.7 days, with amplitudes 
ranging from 1.0-140.8 mmags. The rotation period distribution is clearly bimodal, 
with peaks at ~ 19 and ~ 33 days, hinting at two distinct waves of star formation, 
a hypothesis that is supported by the fact that slower rotators tend to have larger 
proper motions. The two peaks of the rotation period distribution form two distinct 
sequences in period-temperature space, with the period decreasing with increasing 
temperature, reminiscent of the Vaughan-Prcston gap. The period-mass distribution 
of our sample shows no evidence of a transition at the fully convective boundary. On 
the other hand, the slope of the upper envelope of the period-mass relation changes 
sign around 0.55 Mq, below which period rises with decreasing mass. 

Key words: stars: rotation, methods: data analysis, stars: evolution, stars: low-mass, 
stars: magnetic field 



1 INTRODUCTION 

Of the readily observable properties of stars, the rotation 
rate is one which evolves significantly on the main sequence: 
intermediate- and low-mass stars are thought to spin down 
throughout their lifetimes, losing angular momentum via a 
magn etised wind that is linked to their outer convection 
zone (|Kawaler|[l98Sl ; iBouvier et all 1 19971 V Measuring rota- 
tion rates for large numbers of stars over a wide range of 
masses and ages is a long-standing goal in stellar astronomy, 
not only to understand the physical mechanisms driving the 
wind and the resulting angular momentum loss, but also 
to calibrate the relationship between period P I0 t, age t and 
stellar mass M, enabling age estimates to be made for indi- 
vidual stars (|Kawalerlll98l lB arncs 2003). However, period 
measurements for main-sequence stars with well-determined 
ages remain scarce, particularly at low masses, and gyro- 
chronological ages are thus restricted to a li mited range of 
masses and ages and remain very uncertain (|Barnesll2007l l. 
Until the 1990s, stellar rotation rates were mainly mea- 
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sured from spectroscopy, via the rotational broadening of 
absorption lines. These rotational velocity measurements 
provided key insights, particularly the well-known spin - 
down law for Sun-like stars, P TOt oc y/t (|Skumanicblll972h . 
However, these measurements yielded only model-dependent 
constraints on the rotation rate, and were limited to rela- 
tively fast rotators. Using modern wide-field detectors, it is 
possible to measure P Iot directly, by monitoring the bright- 
ness of large numbers of stars simultaneously, and detecting 
quasi-periodic brightness variations, which arise as magneti- 
cally active regions on the star's surface rotate in and out of 
view. Over the past decade, open cluster surveys have pro- 
vided thousands of measurements for low mass stars with 
ag es up to ~ 650 Myr, an overview of which can be found 
in llrwin fc Bouvierl (|2009h . Notable, more recent additions 
to the literature on rotation of early main-sequence low- 
mass stars include studies of M37 (pVleibom et alj|20"09h and 
Coma Berenices l|Collier Cameron et al. 20091 ). Together, 
these data have provided a relatively complete, but complex 
picture of rotational evolution on the pre-main sequence. In 
turn, this has led to renewed efforts to develop models which 
can describe, or even better explain, the observations over 
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the full sub-solar mass range jBarnes fc Kimll20ld ; iBarnej 



l20ld ; lReiners fc Mohantvll2012f ) 



Strictly Periodic 



Correlated Noise 



Period measurements for older stars remain very scarce, 
because they rotate more slowly and are less active than 
their younger counterparts, making it is very difficult to de- 
tect their rotational modulation from the ground. Notable 
exceptions include 71 main sequence F, G and K s tars ob- 
serve d as part of the Mount Wilson HK project (|Barnesl 
I2003T ), 17 27 mid-F to mid -K stars observed by the CoRoT 
satellite (|Affer et all 120121 ') . a nd 41 low mass (0.1-0.3 M B ) 
stars from the MEarthsurvev ijlrwin et al.l201ll ). The Kepler 
space mission (|Borucki et al.ll2010l ) now offers a unique op- 
portunity to measure rotation periods even for slowly rotat- 
ing, moderately active stars, thanks t o its superior precision 
and long baseline. A previous study bv lHarrison et all (|2012h 
measured rotation periods for 265 stars with T c a ^ 5200 K 
and loggr ^ 4.0 dex observed by Kepler for 1-2 quarters 
through the Cycle 1 Guest Observer program. 

The few open clusters included in Kepler's field-of-view 
are particularly important, since their ages can be estimated 
relatively well, and a small sample of perio ds has already 
been published for 71 members of NGC 6811 jMeibom et alj 
l201lt > . Kepler also observed tens of thousands of field stars, 
which can yield period measurements. Although they lack 
individual age estimates, they provide a global picture of 
stellar spin across our Galactic neighbourhood, and can be 
used to constrain the period-mass-age relation in a statis- 
tical sense. In this paper, we focus specifically on the Ke- 
pler M-dwarfs with mass range (0.3-0.55 Mq) where there 
are extremely few previous period determinations for main- 
sequence objects. 

Standard approaches to period detection in light curves 
are based on Fourier decomposition or, for irregularly sam- 
pled data, leas t-squares fitti n g of sinusoidal models and var i- 
ants thereof (|Scargld [l982l ; IZechmeister fc Kiirsterl l2009t) . 
However, typical stellar light curves are neither sinusoidal 
nor strictly periodic, probably because of the clumpy and 
time-evolving nature of the underlying active region distri- 
bution. Residual instrumental systematics are often present 
as well. 

These effects can all lead to a complex periodogram 
structure, with spurious peaks from jumps and long term 
systematics, and multiple or split peaks from spot evolution 
or differential rotation. It is therefore challenging to deter- 
mine which peak corresponds to the rotation period, without 
a priori knowledge of the range of rotation periods expected. 
Consequently, Fourier-domain methods are not always the 
best suited to make the most of Kepler's many thousands of 
spot modulated light curves, which display a wide range of 
rotation periods. We present an alternative approach based 
on the autocorrelation function (ACF) of the light curves. To 
our knowledge, this is the first time that the ACF is used as 
the primary tool t o detect stellar rotation periods, although 
lAffer et all |2012l ) used it as a secondary verification tool. 

In Section [2] we introduce the autocorrelation function 
(ACF) as a robust method for period detection in time series 
data. The results for the Kepler M-dwarf sample are shown 
in Section [31 and discussed in detail in Section 0] Section [S] 
summarises our conclusions and outlines plans for future 
work. 




Figure 1. Simple synthetic signal of amplitude 1 (top row), and 
corresponding ACF (centre row) and periodogram (bottom row). 
The right hand column shows the effect of introducing noise of 
amplitude 0.9, correlated on a 2 time unit timcscale. On the ACF 
and the periodogram panels, the input period used to generate 
the signal is shown as a vertical dotted line, detected period from 
each corresponding method is marked as the over-plotted dashed 
line. 



2 THE ACF METHOD 

In signal processing, the ACF takes the standard form 



E«=i - x) 



(1) 



(see e.g. IShumwav fc Stofferll20icil ) where is the autocor- 
relation coefficient at lag k, for time series Xi (i — 1, . . . , N). 
Each lag k corresponds to r% — kAt, where At is the ca- 
dence. In our implementation, the light curves are median 
normalised before the ACF is computed, and we only search 
for periods less than half the length of the dataset, i.e. k < 
N/2. 

We compare the ACF method to the method most com- 
monly used to search for rotation periods in stellar light 
curves, namely le ast squares fitting of sinusoids over a grid 
of tr ial periods jlrwin et alj 120061 ; IZechmeister fc Kiirsterl 
120091 ). The amplitude, phase and zero-point of the sinusoid 
are free to vary. The sine-fitting periodogram is expressed 
in terms of the statistic 



1 - 
(Xo 



x) /Xo ■ 



(2) 



where Xo is the reduced chi-squared of the light curve with 
respect to a constant value, and \ 2 is the reduced chi- 
squared with respect to the best-fit sinusoid. This m ethod 
is described in more detail in iMcQuillan et all \2Qlj ). 

Fig- [U shows two synthetic time-series curves, together 
with their ACFs and least-squares sine curve fitting peri- 
odograms (|Zechmeister fc Kiirsterl |2009| ) . The left column 
shows a strictly periodic signal, for which the periodogram 
displays a clear pronounced peak. The ACF displays an os- 
cillatory behaviour, with regularly spaced peaks located at 
multiples of the period. The amplitude of these peaks decays 
gradually because of the definite duration of the time-series. 
The right column shows the effect of introducing correlated 
noise to the signal in the left hand column. 

In Fig. [2] we demonstrate how the ACF and the peri- 
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Figure 2. Synthetic examples showing the effects of varying signal phase and amplitude and injecting noise and systcmatics into 
the sine curves (top row), on the ACFs (middle row) and pcriodograms (bottom row). The 'Multiple Effects' signal comprises of the 
'Phase/Amp. Change' signal from the first column, with injected white noise, correlated noise, a linear trend and a jump. On the ACF 
and the periodogram panels, the input period used to generate the signal is shown as a vertical dotted line, detected period from each 
corresponding method is marked as the dashed line, overlaying the dotted line in most cases. 



odogram are affected by varying signal phase and amplitude, 
systematics effects and noise. 

In cases where the phase and amplitude of the signal 
vary with time, the correct period is detected but the ACF 
peak amplitude varies within an envelope, corresponding 
to the amplitude variation of the signal. The peak width 
and level of symmetry can also vary. In this case, the pe- 
riodogram produces two peaks on both sides of the correct 
period. 

If the signal contains multiple minima and maxima per 
period, as can occur for spotted stars with more than one 
dominant active region, the ACF often shows alternating 
low and high ACF peaks, due to a partial correlation be- 
tween the sets of maxima or minima. Our algorithm is built 
to identify these cases (see Section \2.1\ , and in this case se- 
lected the right period. On the other hand, the periodogram 
picked half the right period. 

A jump or long term trend in the signal introduces a 
long term trend in the ACF, and since this can take many 
shapes, it is important to look at the local variations in 
peak height when performing diagnostics on the ACF. These 
long term trends introduce a long-period peak in the peri- 
odogram, which can lead to a wrong identification of long 
periodicity. We therefore expect the ACF method to be more 
reliable in these cases. 

Noise and systematics are present in many stellar light 
curves displaying rotational modulation, and must be ac- 
counted for when attempting to determine rotation peri- 
ods using the ACF. The effect of combining these factors is 
shown in the right column of Fig. [2] which has phase and 
amplitude modulation, white and correlated noise, a linear 
trend and a jump. 

In summary, because the ACF measures only the degree 
of self-similarity of the light curve at a given time lag, the pe- 
riod remains detectable even when the amplitude and phase 



of the photometric modulation evolve significantly during 
the time-span of the observations. The ACF method is also 
capable of producing robust results in cases with residual in- 
strumental systematics, because correlated noise, long-term 
trends and discontinuities give rise to monotonic trends in 
the ACF, on top of which we are able to identify the lo- 
cal maxima. Therefore, the ACF method is expected to 
be more robust to active region evolution than the peri- 
odogram, which implicitly assumes a stable, sinusoidal sig- 
nal. The analysis of the real Kepler data supports this con- 
clusion, as detailed in Section [3JJ] 

2.1 Measuring periods from the ACF 

The period measurement involves three steps: identifying 
peaks in the ACF, selecting the peak associated with the 
mean rotation period, if any, and evaluating the uncertainty 
on the period. 

The presence of high-frequency noise in the light curves 
leads to numerous local extrema in the ACF. Therefore, we 
first smoothed the ACF by convolving it with a Gaussian 
kernel, of window size of 56 lags and a full width at half- 
maximum (FWHM) of 18 lags. These values were tuned to 
provide the best compromise between reducing noise and 
maintaining ACF signal, without prior knowledge of the pe- 
riod. See Fig.[3]for an example of the effects of the smoothing 
treatment. We then identify local extrema in the smoothed 
ACF, defined as locations where the gradient changes sign. 

If the light curve contains a clear rotational modulation 
signal, this process yields a series of clear, regularly spaced 
peaks of gradually decreasing height, as seen in Fig. [TJ The 
first peak corresponds to the interval between patterns in the 
light curve, which evolve gradually, but are clearly repeated, 
and is thus identified as the rotation period. Some of the 
light curves contain long term trends and discontinuities, 
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Tk 



Figure 3. Example of an unsmoothed ACF (black curve) and 
corresponding smoothed curve (red). The many small peaks de- 
tected in the unsmoothed curve are shown as grey vertical dashed 
lines and the peaks detected from the smoothed version are shown 
as red vertical dashed lines. The inset plot shows a section of the 
ACF in detail. 



as a result of imperfections in the systematics correction. 
These introduce power at low frequencies, and thus affect 
the behaviour of the ACF for large lags, but the first ACF 
peak still corresponds to the correct period (as identified by 
visual examination of the light curve). 

In the remaining steps, we make use of the height of 
the ACF peaks. However, correlated noise and residual sys- 
tematics can introduce underlying long term trends, which 
mean the absolute peak height is no longer a good diagnos- 
tic. To mitigate this effect, we measure the height of each 
peak relative to the two adjacent minima, and adopt the 
mean of the two measurements as the 'local height' of the 
peak, denoted by hp. Since the ACF values range from -1 
to 1, hp has a positive value, of maximum 2. 

If a star has two dominant active regions located on 
opposite hemispheres, each causes a series of dips in the light 
curve, approximately in anti-phase with each other. This 
gives rise to a partial correlation at half the mean rotation 
period, leading to a peak in the ACF at Tk = P/2. However, 
the peak at P/2 is typically smaller than that at P, because 
the two active regions do not generally give rise to identical 
dips. We therefore adopt the following condition: if hp of 
the second peak is greater than that of the first, the second 
peak is selected instead. The right panel of Fig. [4] and the 
top left panel of Fig. [5] show a synthetic and real example of 
this effect, where this method has detected the correct peak. 

In some cases, correlated noise and residual systematics 
produce an underlying slope at small Tk , causing a shift the 
position of the first peak associated with the rotation period. 
This occurs because a peak on a slope will have its maxima 
shifted in the direction of increasing gradient, as seen in the 
left panel of Fig. 3] To avoid this bias, a more robust period 
measurement is obtained using the median of the intervals 
Art between consecutive ACF peaks associated with the 
rotation period. To avoid selection of erroneous peaks, only 
those located at or close to (within 20%) of integer multi- 
ples of Tk of the selected peak are used. If there are several 



locations selected around each peak (i.e. the smoothing has 
failed to remove all erroneous peaks) , the peak selection only 
occurs after a gap of At*, > 0.3 Tk- This ensures only 1 data 
point per peak is used in the period measurement. Since the 
accuracy of the peak positions can decrease for very large 
Tk , a maximum of 10 peaks was selected for measurement of 
the median period and uncertainty. 

We define the period uncertainty as the scatter of the 
Arfc. Specifically, we used 

1.483 x MAD 



where N is the number of peaks, and MAD is the median of 
the absolute deviations from the median Ar^. This 'MAD- 
estimated scatter' is equivalent to the standard deviation 
for a Gaussian distribution, but is more robust to outliers. 
In cases where we identified only one peak matching the 
selection criteria, we adopt the peak position as the period 
and the half width at half local peak height as the period 
uncertainty. 



2.2 Tests with Realistic Light Curves with 
Injected Simulated Stellar Modulation 

To demonstrate the robust nature of the ACF, we present a 
selection of synthetic examples with known input parame- 
ters. A more detailed evaluation of the method over the full 
range of parameter space will be described in future work. 

To generate synthetic light curves we used a simple spot 
model code, described in lAigrain et al.l ||2012t ). which takes 
the following input parameters: number of spots, light curve 
amplitude, characteristic spot half-life (A), light curve du- 
ration, fractional differential rotation, inclination and time 
sampling. 

The observations o f lJackson fc Jeffries! (|2012T l show that 
activity is better explained by a large number of small spots, 
than a small number of large spots. They suggest a spot 
number of 2500-5000, however, since they observe younger 
stars we opt for a smaller spot number of 200. This also 
increases the speed of synthetic light curve generation. We 
selected an amplitude of 1% and an inclination angle of 75° . 
The time sampling and duration match that of the Kepler 
Ql-4 data. 

Since the effects of the noise and residual systematics 
of the Kepler data are an important factor when analysing 
the ACF, we introduce these into the synthetic light curves. 
A random selection of the Kepler M-dwarf light curves were 
visually examined and 10 were selected which appeared to 
contain only noise and systematic effects, with negligible 
intrinsic stellar signal. For each synthetic output, we created 
10 light curves with realistic noise and systematics by adding 
each of the selected quiet Kepler light curves. 

The two synthetic examples in Fig. [4] demonstrate the 
ability of the ACF to recover the input period from quite 
difficult cases, where the periodogram failed to identify the 
right period. 

The many synthetic cases we tried led us to believe 
that the ACF method and the periodogram are less reliable 
for light curves with very fast spot evolution (A < 1 period), 
short periods (<7days) and strong systematics. This is 
caused by the steep initial slope in the ACF, combined 
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Figure 4. Examples of synthetic light curves and corresponding ACFs and periodograms for different input parameters, indicated in 
the plot titles. In the left panel the spot distribution was selected to give a double-dip light curve. The dashed lines on the light curve 
plots indicate intervals at the period detected by the ACF. On the ACF and periodogram plots, the input period is marked by a dotted 
line and the detected period is marked by a dashed line. 



with the smoothing algorithm, which masks the peaks. Sim- 
ilarly for fast evolution, strong systematics and long peri- 
ods (> 35 days), the ACF peaks are too strongly affected 
by underlying trends. The periodogram suffers mainly from 
spurious peaks at long periods. 

As the spot half-life is increased to I period, the periodic 
signal strength increases, and by A = 5 periods, the majority 
of periods are correctly detected, across the range of input 
values. It should be noted that these synthetic examples 
were designed to test the limits of the ACF and therefore 
all have a low signal to noise. A recen t study of a sample 
of Sun-like stars observed by CoRoT (|Mosser et all 120091 ) 
found that spot lifetime of the order of the rotation period 
are not atypical for Sun-like stars. On the other hand, the 
spot l ifetime is also known to increase with increased activity 
level (|Hall fc Henrvlll994r i. 



2.3 Amplitude of variability 

To estimate the average amplitude of the periodic signal, we 
use a modified ve r sion o f the range statistic R va , T defi ned by 
iBasri et al.l (|2010l . l201lf ) and lMcQuillan et al.l (|2012f l. which 
is the interval between the 5 th and 95 th percentiles of the 
light curve. This rank-order approach enables us to quantify 
the amplitude of the variability without specifying a par- 
ticular model for it, and excluding the extreme values min- 
imises the sensitivity to outliers and high-frequency noise. 
We divide the light curve into segments, whose duration is 
equal to the detected period, measure R va , T separately in 
each of them, and define R peT as the average of the i? V ar 
measurements obtained in this way. In most cases, -R per is 
very similar to the R va _ r value for the whole light curve. 



3 APPLICATION TO KEPLER DATA 
M-DWARFS 

3.1 Target Selection 

Following ICiardi et al.l (|201ll) . we selected M-dwarfs among 
the Kepler targets based on the effective temperature T e g 
and surface gravity log g report ed in the Kepler Input Cata- 
log, or KIC (|Brown et al.ll20rif ). These parameters were es- 
timated using Baye sian posterior probability maximisation 
(|Brown et al.l 120111 ). matching observed colours, estimated 
from Sloan g, r, i, z filters, 2MASS J HK, and D51 (510 
nm), t o the stellar atmosphere models of ICastelli fc Kurucj 
|2004h . resulting in typical uncertainties of 200 K for T e s 
and 0.5 dex for logg, respectively. These are the typical er- 
rors for each parameter, but the actual errors may vary by 
a small amount between stars of different magnitude and 
spectral type. F or a more detailed discussion of the KIC 
par ameters, see iBrown et all (|201lh . iBatalha et ail (|2010T l 
and IVerner et alT i 201lf ). There is increasing evidence that 
the KIC values for the radii and r e ff are overestimated (e.; 
iMann et al]|2012l ; IVerner et al1l201ll ; iMuirhead et al.ll2oW 
which means the values of T e ff should be interpreted with 
care. 

This study focusses on M-dwa rfs stars, which were se- 
lected following ICiardi et al.l |201lJ ) as having T c ff ^ 4000 K 
and logg ^ 4.0 de x, manually adding 25 known M-dwarfs 
(|Ciardi et al.ll201ll ) with missing KIC parameters. Table Q] 
summarises the number of stars considered at each stage of 
the study. 



3.2 Obtaining and pre-processing the light curves 

This study is based on public release 14, which was avail- 
able when the present analysis was performed. The release 
included data from Quarter 1 (Ql) to Quarter 4 (Q4) of 
Kepler observations, which took place over ~ 310 days be- 
tween May 13 th 2009 and March 19 th 2010. Approximately 
156 000 targets were observed during this time, most with 
a cadence of 29.42 minutes. These data are publicly avail- 
able and were downloaded from the Kepler mission archive 
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Table 1. Number of objects included at each stage of the study. 
The column labelled 'All dM' includes both M-dwarfs selected 
from the KIC parameters and the previously known M-dwarfs 
without KIC parameters. The totals in the second column include 
the known EBs and planet candidate host stars, which are also 
listed separately in the last two columns. The possible giants and 
double-period stars were removed from the final results reported 
in Tables [2] & [5] Note: removed 4 contact EBs from Ql-4 EB 
selection. 



Stage 


All dM 


EBs 


PI. host 


KIC selection 


2937 


15 


57 


Ql-4 light curves 


2483 


9 


51 


Period detected 


1730 


7 


42 


Possible giants 


121 








Double-period 


39 








Binaries or pulsators 


112 







Final periodic total 


1570 


7 


42 



(http://archive.stsci.edu/kepler) at the Space Tele- 
scope Science Institute (STScI). 

Public release 14 included a re-processing of Q1-Q4 
with a new version of the presearch data conditioning (PDC) 
pipeline known as PDCMAP. The purpose of the PDC is to 
remove the majority of instrumental glitches and systematic 
trends. In contrast with earlier versions, the new PDC-MAP 
uses a Bayesian approach t o do so while retaining most real 
(astr ophysical) variability l|Smith et al.ll20l"3 ; IStumpe et alj 
120121 ). It works very well for Q3 and Q4, but some artefacts 
remain in some light curves in Ql and Q2. Improving the 
correction of instrumental effects further remains an impor- 
tant goal to enable the full exploitation of Kepler's potential 
in terms of stellar astrophysics. However, in the vast major- 
ity of the cases we have examined, the residual artefacts in 
the PDC-MAP data do not prevent the detection of rota- 
tional modulation, so we consider it suitable for the present 
work. 

An additional 5 quarters of PDC-MAP data were made 
publicly available a few weeks before the submission of this 
manuscript. We are planning to analyse these in the near 
future, in the hope that it will improve our sensitivity to 
long periods, but we do not expect it to affect the main 
conclusions of this paper. 

The ACF calculation requires the light curves to be reg- 
ularly sampled and normalised to zero. We divided the flux 
in each quarter by its median and subtracted unity. Gaps 
in the light curve longer than the Kepler long cadence were 
filled using linear interpolation with added white Gaussian 
noise. This noise level was estimated using the variance of 
the residuals following subtraction of a smoothed version of 
the flux. To smooth the flux we applied an iterative non- 
linear filter which consists of a median filter followed by a 
boxcar filter, both with 11-point windows, with iterative 3- 
sigma clipping of outliers. 

This method of gap filling and quarter stitching was 
only found to introduce spurious jumps in light curves for 
which the systematics correction had failed, and were there- 
fore already problematic. In the vast majority of cases, the 
amplitude of M-dwarf activity mean that any quarter stitch- 
ing effects have a negligible impact on the period detection. 

Our sample also includes a number of previously iden- 



tified eclipsing b i naries and planetar y transit candidates 
(|Pr§a et all |201ll : iBatalha et all l2012t ). In those cases we 
cut out and interpolated over the eclipses/transits before 
performing the rotation period search. 

3.3 Establishing confidence in the period 
detection 

Given the scientific importance and manageable size of the 
sample under study here, we opted to perform a visual exam- 
ination of all the light curves in order to verify the period 
detected by the ACF. The light curves were compared to 
dashed vertical lines at intervals of the detected period (see 
e.g. Fig. [5]l . In order to classify the detection as valid, fea- 
tures in the light curve must be present at intervals matching 
the dashed lines, across several periods, preferably in more 
than 1 quarter. 

We paid particular attention to the question of whether 
the detected period was clearly the rotation period, or could 
be P/2, as described in Section 12.11 In 73% of light curves 
we identified features which could be 'tracked' visually, i.e. 
a particular shape of spot crossing that repeats throughout 
the light curve. When there are more than one such sets of 
features, and they evolve gradually relative to each other, as 
shown for example in the top left panel of Fig.0 the rotation 
period is very clear. When only one such set is visible and 
the first ACF peak is not higher than the second (e.g. the 
bottom right panel of Fig. [SJ, we cannot be so certain, but 
the simplest explanation is nonetheless that the detected 
period is correct. In around 1% of the cases where we report 
a period, there is very tentative evidence in the light curve 
that the detected period is a harmonic (i.e. P/2), but we 
could not be certain, so we simply flagged the corresponding 
objects as 'possible harmonics'. Such cases may be resolved 
once additional quarters of Kepler data are incorporated into 
the analysis. We note that these cases are not numerous 
enough to explain the two sequences seen in the period- 
temperature diagram (Fig. 1101) . and are spread across the 
entire period range. 

In 3.6% (57) of the periodic cases an incorrect peak 
was identified. This was most frequently a result of noise 
introducing extra peaks, or very large residual systematics 
changing the relative peak heights. These cases were man- 
ually corrected to select peaks corresponding to the period 
identified by eye. 

In the future, we plan to apply the same analysis to a 
much larger sample of Kepler targets (including F, G and 
K dwarfs), for which an automated detection method must 
be developed. The present, visually inspected subset will 
then prove valuable, to validate any threshold to be applied 
automatically to a larger sample. Here we merely note that, 
in the present study using 10 months of data, it is possible 
to recover 91% of our detections, at the cost of a false alarm 
rate of 10%, by selecting objects with 

hp > MAX (°- 15 ' 51^) • (4) 

3.4 Results 

We applied the ACF method and visual inspection 
steps described in the previous section to the Ke- 
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Figure 5. Examples of Kepler light curves and corresponding ACFs. The dashed lines on the light curve plots indicate intervals at the 
period detected by the ACF (dashed line on ACF plot). KIC 4918333 shows two active regions, creating a double dip effect. The correct 
period was automatically detected by selecting the second (higher) ACF peak. KIC 9201463 shows that the ACF is robust against flares 
and significant systematics. KIC 9488191 contains significant residual systematics and the effect on the ACF is apparent, however, the 
correct period is still detected. KIC 7590933 shows an example which has been marked as 'possible harmonic' or 'HM' in Table [2] since 
it is not clear whether the rotation period corresponds to the first or second ACF peak, even though the ACF peak height indicates it 
should be the first (see text for further discussion). 



Table 2. M-dwarfs with detected rotation periods. This table is available in its entirety, in a machine-readable form in the online 
supplementary material, or at http://www.physics.ox.ac.uk/StellarRotation A portion is shown here for g uidance regarding it s form 
and content. T e g and logg are from the KIC and M was derived from T e ff using the 600 Myr isochrone of iBaraffe et alj dl998f) . The 
average amplitude of var iability per period bin of the li ght curve, R ver , is in cluded. The meaning of the flags are: 'EES': known eclipsing 
binary jPrsa et al,ll201ll) ; 'PL': planet host candidate l lBatalha et al J [20121) : 'PB': ultra-stable periodic behaviour, indicating possible 
binary, pulsator or young object; 'HM': the reported period may be a harmonic of the true period; 'NF': no flag. 
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0.046 


10.4 


NF 


1572802 


3990 


4.48 
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Table 3. Objects identified as likely giants from their light 
curves. This table is available in its entirety, in a machine- 
readable form in the online supplementary material, or at 
http://www.physics.ox.ac.uk/StellarRotation A portion is 
shown here for guidance regarding its form and content. 



KIC 
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log g 








(M ) 


(g/cm 3 ) 


(K) 


(mmag) 


1026895 


0.5343 


4.53 


3977 


8.99 


1160867 


0.4474 


4.57 


3753 


8.29 


1431599 


0.5068 


4.5 


3907 


12.68 


1576043 


0.4889 


4.52 


3861 


7.76 


2010137 


0.4881 


4.48 


3859 


7.01 



pier light curves of the objects selected as likely M- 
dwarfs. The results are reported in Tables [5] to [5] The 
full machine-readable tables are available in the online 
supplementary material, or with plots of every light 
curve, its ACF, and sine-fitting periodogram at the URL 
http : //www .physics . ox. ac .uk/StellarRotation . Table [2] 
reports all our period measurements, except for two groups 
of objects, listed separately in Tables and 2] which were 
excluded from the sample for reasons detailed below. Table[5] 
lists all the objects which passed our target selection criteria, 
but for which no period was detected. The final number of 
likely M-dwarfs with detected rotation periods is 1570 (out 
of the 2483 light curves we analysed, or 2362 if we discount 
objects later found to be giants). 



3.5 Excluding non- rotators 

Rotational modulation of star spots is not the only cause of 
periodic, or quasi-periodic, variability in stellar light curves. 
The visual examination stage was therefore important in 
identifying groups of stars whose variability did not appear 
to be caused by rotation, or by rotation alone. 

The Kepler light curves have al ready been searched 
extensively for pl anetary transits (Borucki c t al.l 20 111 ; 
iBatalha etall 120121 ) and stellar eclipses (|Prsa et al.ll201lf ). 
We included these objects in our sample, after removing 
the transits and eclipses, as described earlier. Measuring the 
photometric rotation periods of known binaries is useful to 
check for differences between the rotational properties of 
close binaries and of apparently single stars. The rotation of 
planet-host stars is an important topic in itself, which is be- 
yond the scope of the present paper and will be investigated 
in detail in a forthcoming publication. 

One distinct group of 121 stars show a clear ACF peak 
at relatively short Tfe (~ 1-12 days), but there are few ad- 
ditional peaks at integer multiples of the first, and the light 
curve appears stochastic rather than truly periodic. Such an 
example is shown in Fig.[6] We initially thought that these 
objects may be rapid rotators with very rapidly evolving 
active regions, or a hitherto unidentified type of pulsating 
star. We examined their KIC parameters, and noted that 
they appear redder, and have lower proper motions, than 
the rest of the sample. We therefore concluded that they 
are likely to be giants contaminating our samp le. Indeed, if 
we ap ply the J — H < 0.75 cut advocated bv lCiardi et al.l 
(|201ll ), it removes 117 stars from the M-dwarf sample, of 
which 103 belong to this group of objects displaying stochas- 




30 40 50 
T k (days) 



Figure 6. KIC 2010137 is an example of a star which we identified 
as a likely giant based on its light curve, which shows stochastic 
variability with a clear, dominant time-scale, but no repeatability, 
and does not resemble the signal expected from rotation. 




Figure 7. Examples of the objects classified as 'PES'. KIC 5516671 
(top) shows a very stable light curve and KIC 3103752 (bottom) 
shows the 'beat pattern' discussed in the text. 



tic behaviour. We therefore concluded that these objects are 
likely giants, removed them from the M-dwarf sample and 
list them separately in Table [3] 

A further 39 objects display evidence for two, very dis- 
tinct periods in their light curves, with one period several 
times longer than the other. This could potentially result 
from the light of two different variable stars being included 
in a single photometric aperture, which may be revealed by 
examination of the pixel level data. Since the nature of the 
periodicity is unknown, and the determination of the KIC 
parameters used to select these objects as M-dwarfs could 
be affected by the presence of a close companion, we exclude 
these objects from our sample, and list them separately in 
Table[4] Note we do not report periods for them, because the 
presence of two distinct signals in their light curves makes 
the identification of either period more challenging. 

Finally, 109 stars show unusually stable periodic be- 
haviour, with periods typically < 7 days, and very little or 
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Table 4. Objects with two distinct periods in their light 
curves. This table is available in its entirety, in a machine- 
readable form in the online supplementary material, or at 
http://www.physics.ox.ac.uk/StellarRotation A portion is 
shown here for guidance regarding its form and content. 
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Table 5. Objects with no rotation period detection. 
This table is available in its entirety, in a machine- 
readable form in the online supplementary material, or at 
http://www.physics.ox.ac.uk/StellarRotation A portion is 
shown here for guidance regarding its form and content. 
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no evidence of any evolution over the full Q1-Q4 duration. 
Some of these also show a beat pattern characteristic of the 
light curves containing two, mutually similar periodicities. 
See Fig. [7] for examples. These objects could be members 
of close binary systems, where the active region pattern is 
stabilised over long timescales because of the presence of a 
companion. Their rotation periods and amplitudes are cer- 
tainly similar to those of the known eclipsing binaries in the 
sample. If they are binaries, the beat patterns may indicate 
differential rotation on one of the stars, or two slightly dif- 
ferent rotation periods for the two components of the binary. 
We are not aware of any type of main-sequence M-dwarf that 
would be expected to pulsate in this period range, but we 
cannot rule out binarity or pulsation without spectroscopy. 
We therefore kept these objects in Table(2] but flagged them 
as 'binaries or pulsators', indicating their period may result 
from a phenomenon other than rotation. 



3.6 Comparison of the AFC and sine-fitting 
periodogram 

We computed sine-fitting periodograms for the sample of 
Kepler M-dwarfs using 1000 logarithmically-spaced periods 
between 0.1 and 155 days. 

In general, similar results are obtained from the ACF or 
periodogram methods, however, we consider the ACF clearer 
and more reliable. A useful feature of the ACF method is 
that it enables automatic identification of cases where the 
detected period is a half the mean rotation period. An exam- 
ple of this is shown for KIC 4918333 in Fig. O which clearly 
shows that the mean rotation period corresponds to the sec- 
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Figure 8. Comparison of the periods detected using the ACF 
and the periodogram, for all stars with a visually confirmed pe- 
riodicity and ACF detection. The correct ACF detections where 
the periodogram detections are within 10% of the ACF period 
are shown in black (1298), the cases where the periodogram pe- 
riod is not within 10% of the ACF period are shown in dark 
blue (215). Cases where the ACF initially selected the wrong pe- 
riod but the periodogram selected the correct one are marked 
by red points (45). Cases where both the ACF and periodogram 
detected the wrong period are marked as cyan points (12). The 
grey dashed lines show where Pacf = PlS; Pacf = 2 x Pls an d 
Pacf = 0.5 x P LS . 



ond ACF peak, whereas the first periodogram peak is the 
highest. 

Similarly, for KIC 9488191 in Fig. [5] one can determine 
easily from the ACF that the rotation period corresponds 
to the first peak at ~ 15 days and not the highest peak in 
the periodogram at ~ 70 days. Although long period peaks 
could be excluded from the periodogram to allow the correct 
peak to be selected, this would require prior knowledge of 
the possible range of rotation periods to avoid removal of 
genuine long period signatures. 

A more systematic comparison between the ACF and 
periodogram results is shown in Fig. [8] Here we make a one- 
to-one comparison of the periods detected using the ACF 
and the sine-fitting periodogram, for all the cases where 
the light curves were determined to show periodic variations 
from visual examination. 

In 222 (14%) of the periodic cases, the periodogram 
does not detect the correct period (to within 10%). The 
most common discrepancy arises from cases where the pe- 
riodogram has detected half the period of the ACF (170 
cases). We stress that there is no obvious way of identify- 
ing these cases automatically using the periodogram alone 
(for example, we tried to use the relative heights of the 
periodogram peaks, without success). In 11 cases the pe- 
riodogram selects erroneous long period peaks and the ACF 
selects the correct shorter period peak. 

On the other hand, the ACF yielded 57 incorrect ACF 
period detections (3.6%), 45 cases occur where the ACF 
method selects the wrong period, and the periodogram se- 
lects the correct period. These cases are marked as red points 
in the Figure [8] For the remaining 12, both the ACF and 
periodogram detect the wrong period (cyan points). 
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Figure 9. Period versus amplitude for the rotating Kepler field M-dwarfs. Blue dots represent objects with P ro t < 1 days, which als o 
display unusually stable modulation patterns in their light curves, a nd blue stars known, short-period eclipsing binaries l lPrsa et al.l20lfl) . 
Red dots represent the host stars of candidate transiting planets dBatalha et al.lfeoi^ . All the other M-dwarfs with detected rotation 
periods are shown as grey dots. The histograms of each parameter are shown along the corresponding axis, with matching colours. Two 
long-period binaries are not shown as blue stars in the plot. 
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Figure 10. Period versus effective temperature for the rotating 
Kepler field M-dwarfs. The symbols and colours as the same as in 
Fig. |5] The apparent dearth of objects around T e ff = 3775 K is 
not apparent in any of the KIC colours, and is therefore thought 
to be a result of the KIC stellar parameter estimation procedure, 
rather than a real effect. The dashed line marks the location of 
the cut made by eye between the fast and slow rotators. 

In one case (KIC 10553513, not included in the periodic 
sample), the large residual systematics prevent either the 
ACF or the periodogram from detecting the rotation period 
which is visible in the light curve. The only case for which the 
periodogram is able to detect the period and the ACF is not 



(even with manual correction) is KIC 5480340 (not included 
in the periodic sample). This is due to the extremely short 
period (0.25 days). The steep gradient in this region of the 
ACF prevents peak detection, which we plan to resolve in a 
future version of the algorithm. 

We conclude that the periodogram remains a valuable 
technique for period detection, however the clarity and ro- 
bustness of the ACF method makes it our tool of choice for 
the measurement of stellar rotation in light curves. The ACF 
method works independently of periodogram based meth- 
ods and is ideal for determining rotation period statistics 
for large datasets. 



4 DISCUSSION 

4.1 Period, amplitude and temperature 

We now examine the period distribution of our sample as 
a function of amplitude (Fig. O and effective temperature 
(Fig. [10]). A number of interesting features are immediately 
apparent. 

First, the period distribution is clearly bimodal for 
Prot > 10 days, with peaks of approximately equal height 
at ~ 19 and ~ 33 days. This bimodality appears statisti- 
cally significant in log period space, with a Hart igan's dip 
test p- value of 0.0003 (|Hartigan fc Hartiganll 19851 ). The two 
peaks of the rotation period distribution form two distinct 
sequences in period-temperature space (Fig. 1 10|) . The pe- 
riod decreases with increasing temperature, in much the 
same way, for both sequences. This variation as a function 
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Figure 11. Histogram of the short and long period M-dwarfs, as 
defined by the cut marked with the dashed line in Fig. 1101 The 
dashed histogram marks the long periods divided by 2, and pro- 
vides further evidence that the short period set are unlikely to 
be erroneous half-period measurements of the long period popu- 
lation. 



Figure 13. Histogram of the total proper motion for the short 
period and long period M-dwarfs, as defined by the cut marked 
with the dashed line in Fig. 1101 The uncertainty on the proper 
motion values is 0.02" /year but given the high numbers in each 
sample these results are still significant. The known and 'possible' 
binaries have been removed from this sample (blue points and 
stars in Fig. 1101 . 
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Figure 12. Total proper motion against Period for the short 
period (green) and long period (black) M-dwarfs, as defined by 
the cut marked with the dashed line in Fig. 1101 



of temperature implies the gap between the sequences is not 
a result of systematics in the Kepler light curves leading to 
missed detections in a particular period range. We define 
the short and long period samples based on the line, plotted 
somewhat arbitrarily, shown in Fig. [TUJ 

Before interpreting this result, we must address the pos- 
sibility that the bimodal distribution is spurious: an error 
of a factor of two in the periods of about half the objects 
could give rise to the observed distribution. Such errors are 
not uncommo n in rotational studies based on ground-based 
data (see e.g. ICollier Cameron et al.l [20091 ). However, they 
are less likely to be prevalent in the present study, for the 
following reasons. We can exclude a scenario where we mea- 
sured twice the true period for the objects belonging to the 
longer period peak, because of the continuous sampling of 
our data. The alternative scenario is that we measured half 



the true period for the objects belonging to the shorter pe- 
riod peak. This kind of problem arises when the brightness 
distribution of the stellar surface is bimodal, due to concen- 
trations of active regions on opposite hemispheres, resulting 
in 'double-dip' light curves. However, our ACF peak selec- 
tion routine was specifically designed to address this prob- 
lem - this is one of the strengths of the method. In fact, for 
210 of the objects belonging to the shorter period peak, a 
'double-dip' light curve lead to the second ACF peak being 
automatically selected as the rotation period. Unless these 
stars have quadrupolar surface brightnesses, it is very un- 
likely that we underestimated their periods by a factor of 
two. 

The case for genuine bimodality is further strengthened 
by comparing the detected distribution to that which would 
arise if either set were wrong by a factor of 2. Fig. 1111 shows 
that division of the long period set by 2 does not reproduce 
the short period set (and equivalently neither does multipli- 
cation of the short period group). A quantitative two-sample 
Kolmogorov-Smirnov (KS) test of the short period sample 
and the long period sample divided by 2 gives a p-value 
9 x 10 -14 , confirming they are very unlikely to be drawn 
from the same distribution. 

Fig. U2l shows a weak correlation between period and to- 
tal proper motion for the two samples. We also see tentative 
evidence of a difference in the proper motion distributions 
of the short and long period samples. Fig. 1 131 shows the his- 
tograms for stars in the short and long period samples with 
non-zero proper motion KIC measurements (266 short, 345 
long). The proper motion values from the KIC are taken 
from a selection of catalog^] where available. Total proper 
motion is listed on NStED as having accuracy of 20 milliarc- 



1 Kepler Stellar Classification Program, Hipparcos, Tycho-2, 
UCAC2, 2MASS and USNO-B1.0. 
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seconds per year, but the large number of objects in each 
sample support the difference between them. 

We also checked for differences in the galactic latitude 
and Kepler magnitude distributions of the two samples but 
did not see any clear differences. By using the galactic lat- 
itude measurements and an approximate distance estimate 
from average apparent and observed M-dwarf magnitudes, 
we conclude that the observed population lies within a small 
fraction of the scale height of the thin disk, and therefore 
a variation of properties with Galactic latitude is not ex- 
pected. 

For the remainder of this study, we therefore assume 
our sample contains two distinct stellar populations. The 
gap between the two sequences in Fig. \W\ is reminiscent 
of the Vaughan-Preston ga p (hereafter, V-P gap) found by 
l|Vaughan fc Prestonlll98Cl . hereafter, VP80) when studying 
the chromospheric activity levels of a relatively small sample 
of F and G stars as a function of B — V colour. As chromo- 
spheric activity, like stellar rotation, declines with age, both 
observed gaps suggest two waves of star formation in the 
Solar neighbourhood, as discussed by VP80 themselv es and 
by numerous authors subseq uently (see, for example, iBarrvl 
Il988l . [Pace fc Pasquinil [20041 and references therein). Alter- 
native explanations for the V-P gap include a discontinuity 
in the chromospheric activity level or of the rotation pe- 
riod a t a particular point in the star's evolution. iNoves et al.l 
found no evidence for a difference in the relationship 
between chromospheric activity, rotation period and colour 
for stars on either side of the V-P gap, which led them to dis- 
miss the first of these alternatives. Our results also suggest 
the gap does not result from a discontinuity in the chromo- 
spheric activity level, as we see no obvious difference in the 
amplitude of photometric variability of the slow and fast ro- 
tators. Of the two remaining hypotheses, that of two waves 
of star formation seems more plausible, and is supported 
by the tentative differences we observe between the proper 
motion distributions of the two groups, but more kinematic 
data would be needed to confirm different epochs of star 
formation. 

The fact that these have different median periods and 
different proper motion characteristics suggests that they 
may result from two distinct waves of star formation: as al- 
ready discussed, low-mass main-sequence stars spin down 
as they age, and older populations tend to have larger ve- 
locity dispersions than younger populations, due to dynam- 
ical heating of the Galactic disk over long timescales (see 
iFreeman fc Bland-Hawthornll2002l . and references therein). 
The ratio of the two peaks in the period distribution corre- 
sponds to an age ratio of ~ 3, if one assumes th at M-dwarfs 
spin down as i 0,5 , as observed for Sun-like stars (|Skumanichl 
Il972h . Recent rotational studies of low-mass stars in open- 
clusters suggest that the index of the spin-d own low might 
be closer to 0.6 (see e.g. iMeibom" et al.ll2009l ). which would 
imply an age ratio of ~ 2.5. We note that an interpretation in 
terms of thin and thick disk populations is unlikely because 
there should be many fewer thick disk than thick objects in 
the Kepler sample, which is essentially magnitude-limited. 

Further kinematic and distance information is required 
to draw stronger conclusions about the nature of the two 
populations. The Kepler mission itself may, in the medium 
term, provide this information: the pixel-level data can be 
used to monitor the centroid of each star over the full 



lifetime of the mission, which should enable the measure- 
ment of pro p er m otions and parallaxes for some targets. 
iMonet et all (|2010t ) shows that astrometric precision for a 
single 30 minute measure is < 4 milliarcsecond s. In the 
longer term, the GAIA mission (|de Bruiineil20f 2l ) will pro- 
vide this information with greater precision and for more 
objects. 

While the bulk of the objects (89%) have periods in 
the range 10-50 days, a small fraction (8%) of the periods 
are < 7 days. Almost all these rapid rotators display unusu- 
ally stable light curves over the full 10 month dataset (blue 
points in Figs. [9] and [10] label 'PB' in Table[2}. Most of them 
also have relatively large amplitudes (> 2%, compared to 
0.2-2% for the bulk of the sample). Two possible explana- 
tions spring to mind for these rapid rotators: they may be 
significantly younger than the rest of the sample (for exam- 
ple members of the young disk population, see Section 14.21) . 
or they may be close binary systems, which have become 
rotators due to spin-orbit interactions. The 5 known short- 
period eclipsing binaries in our sample (shown as stars in 
Fig. [9} all display synchronised rotation signals, suggesting 
that at least some of the other However, we cannot con- 
clusively distinguish between the two possibilities without 
spectroscopy. 

Finally, all the planet-host candidates (jBatalha et al.l 

|2012| , red points in Figs. [9] and [10} have P rot > 10 days and 
amplitude < 2%. If the period and amplitude distributions 
of the planet-host candidates were the same as those of the 
rest of the sample, we would have expected ~ 3 of them to 
have P r ot < 10 days and/or amplitude > 2% (of all the M- 
dwarfs with rotation periods, 10% fall outside these limits, 
and there are 42 candidate planet host stars in the sample) . 
The number of objects concerned is too small to draw firm 
conclusions, but it does suggest that either the search for 
transits in Kepler light curves may be less complete around 
active, rapidly rotating stars, or that there are fewer tran- 
siting planets around these stars (which could be the case 
of the rapid rotators are close binaries). Either way, this 
should be taken into account when inferring the incidence 
of planets around different types of stars. 

4.2 Period-mass relation 

To compare our results with other rotational studies of low- 
mass stars, we estimated masses for the Kepler M-dwarfs 
from the KIC effective temper atures. Much of the previous 
work on rotational evolution (IKaw alcr 19 891: iBarnesI 20031 . 
120071 : IMeibom et"al"1l2009l : ICollier Cameron et al.ll2009l ) uses 
a directly observed colour index, such as B — V or J — K, in- 
stead of model-dependent mass estimates, to confront mod- 
els with observations. We opted for masses for three rea- 
sons. First, angular momentum evolution models generally 
depend on fundamental parameters rather than colours. Sec- 
ond, neither of the colour indices most frequently used, B—V 
and J — K , is well matched to our targets, which are very 
faint in B, and have almost constant J — K over the mass 
range of interest (in contrast to G and K-dwarfs, where J—K 
is a steeper function of mass). Third, stellar parameters 
based on multi-colour photometry should be less sensitive 
to reddening than any single colour index. 

The masses used in this stud y were obtained by inter- 
polating the 630 Myr isochrone of iBaraffe et al.l (|l99Sl ). We 
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Figure 14. Period ver sus mass for low-m ass field stars. Based on data from lBaliunas et alj l|l996h and lKiraga fc St cpicn (20q3) (circles) 
and MEarth data from llrwin et all j201ll ) (stars), with grey and black symb ols representing object s with young and old disk kinematics, 
respectively. Additional M-dwarf periods from the WFCAM Transit Survey llGoulding et alj|2012h , for which no kinematic classification 
is available, are shown as white triangles. The new results from the present study are shown as blue points. 



checked that the results are essentially unchanged if the age 
is increased by a factor of up to 10. While useful for a general 
comparison, we stress that these masses are very uncertain 
and should not be used for individual objects. The 200 K 
formal uncertainties on T e g alone translate into mass uncer- 
tainties ranging from 0.15 M Q (for M ~ 0.3 M©) to 0.1 M© 
(for M ~ 0.6Mq). Additionally, any systematic errors in 
T c ff, discussed in Section [3.11 are automatically propagated 
through to the masses. 

Fig. [TJ shows the period-mass relation for field stars 
with masses b elow 1.3 Mb, based on a compilation of liter- 
ature sources llBa iunas et al.l fl996: Ki raga fc Stepienll2007l : 



llrwin et alj|201ll ; iGoulding et al]|2012l ). Our sample (small 
blue dots on Fig. I14|) ties in smoothly with the existing data 
at both higher and lower masses (circles, stars and trian- 
gles). The tight sequence followed by old F, G and K stars 
is clearly seen on Fig. 1141 as are the younger stars gradually 
evolving towards that sequence. This evolution is relatively 
well reproduced by sim ple models of angular momentum loss 
via a magnetised wind l|KawaIeilll98i , ll989h . and forms the 
basis for ro tational age-dating, or gyrochronology (|Barnes! 
I2003l . l2007l ). 

The picture is more complex for M-dwarfs, which evolve 
more slowly, and hence have not converged onto a common 
sequence, even after several Gyr. However, the new, much 
larger sample of early M-dwarfs strengthens two interesting 
features, which were only hinted at by previously published 
results. 

First, we might h ave expected to see some kind of transi- 
tion a round 0.35 M© (|Chabrier fc B araffc 1997; S cholz et alj 
l201ll ). which corresponds to the transition between fully con- 
vective stars and stars with a radiative core, because the so- 
called interface dynamo, which powers the magnetic field of 
Sun-like stars, cannot operate in the absence of a radiative 
core. No such transition is seen among the slowly rotating 



stars, which define the upper envelope of the period-mass 
relation. There does seem to be a transition at ~ 0.35 M© 
for rapid rotators: below this mass, rapid rotators persist 
e ven among k i nema tically old stars, as the MEarth sample 
of|l rwin et all (|201 lh illustrates, b ut above it, all the f ast ro - 
tators are kinematically young, as lKiraga fc Stepienl (|2007l ) 
pointed out. The Kepler sample as it stands cannot shed 
much additional light on this point because we currently 
lack kinematic and distance information. 



Second, a transition that was not expected, but is 
clearly seen, occurs in the upper envelope of the period-mass 
relation around 0.6-0.55 Mq, where the slope of the relation 
suddenly changes sign. To our knowledge, this intriguing 
feature has not yet been discussed in the literature, and it 
will need to be accounted for in future modelling work. In 
particular, this transition will need to be incorporated into 
empirical gyrochronological relations, if the latter are to be 
applied to M-dwarfs. 



The moderately slow rotation of field M-dwarfs has been 
problematic for some time. The latest generation of theo- 
retical moA^_^_mi^ l \A^_mor^nt\mi loss via a magnetised 
wind (|Reiners fc Mohantv||2012l ) naturally explains the wide 
range of rotation rates observed for M-dwarfs in early main- 
sequence clusters. However, the behaviour of field M-dwarfs 
can only be reproduced by assuming that the critical rota- 
tion rate (above which the magnetic field saturates) is not 
universal, but depends on mass and perhaps even on age. 
With the new results, upper-envelope of the period-mass 
relation is now much better defined, which should prove 
valuable in testing and calibrating future refinements of the 
models. 
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5 CONCLUSIONS 

We report rotation periods for 1570 main-sequence M-dwarf 
stars with masses between 0.3 and 0.55 Mq, measured using 
a new method based on the autocorrelation function, applied 
to the first 10 months of data from Kepler. The fraction of 
objects in which we de tected periods, 63 .2%, is remarkably 
high. For comparison, llrwin et al.l (|201lh detected periods 
in 15% of their 273-object sample. The contrast can be ex- 
plained by the unprecedented precision and continuous sam- 
pling of the Kepler data, and by the performance of the ACF 
method. 

The bulk of our period detections fall into two distinct 
groups, with periods in the range 10-25 and 25-80 days with 
peaks at ~ 19 and ~ 33 days respectively. We suggest that 
these correspond to two stellar populations with different 
median ages. The comparison of available non-zero proper 
motions for the two main samples further supports the hy- 
pothesis that they have different ages. The more slowly ro- 
tating group has typically higher proper motion values, as 
would be expected for an older population. Within each 
group there is also a weak correlation between period and 
proper motion. This study shows the potential use of vari- 
ability statistics as a probe of star formation history, al- 
though further data on distance and kinematics are required 
to draw sound conclusions about the nature of these regions. 

In each of the two groups of stars, which form the two 
peaks of the period distribution, the rotation period tends 
to increase with decreasing mass. This trend is very clear 
for the upper envelope of the relation, but we also checked 
that it applies to the bulk of the stars by examining the 
median of the rotation periods for each group (slow and fast 
rotators) in 0.05 Mq bins. By combining our results existing 
rotation data from the literature over the mass range 0.1- 
1.3 Mq, we have shown that this relation extends over the 
whole M-dwarf regime (0.1-0.55 Mq) but is in stark contrast 
to the behaviour of K-stars, whose period decreases with 
decreasing mass. To our knowledge, this dichotomy between 
K and M-dwarfs has not been noticed before, and is not 
predicted by current models. 

A small fraction (~ 8%) of the M-dwarfs display short 
(< 7days), stable periods, and marginally enhanced vari- 
ability compared to the rest of the sample. The most likely 
explanations are that these are short-period binaries or very 
young stars. 

When combining the Kepler sample with data from the 
literature, we see no evidence for a break in the period-mass 
relation around 0.35 Mq, even though stars below this mass 
are expected to remain fully convective. Such a break might 
have been expected if the development of a radiative core 
played a key role in drivi ng a large-scal e magnetic field, as 
proposed for example by iBarne 1 (|2003h . We note that the 
Kepler and MEarth field M -dwarfs lie in the 'rotational gap' 
defined by iBarnej ^OIOD . This suggests that the Rossby 
number, which plays a key role in controlling the rotational 
evolution of G and K stars, may be less important for low 
mass stars. 

Analysis of additional quarters of Kepler data may re- 
veal long periods for some of the objects without detections 
in the present sample. We are also working to improve the 
systematics correction with an alternative to the PDC-MAP 
(Roberts et al. submitted) , and to enhance the performance 



of the ACF method. Residual systematics and quarter joins 
are the limiting factors in the current analysis because they 
introduce steep variations in the ACF. By optimising the 
ACF smoothing parameters after initial period detection, 
and working on methods to separate the periodic signal from 
the long term ACF trends, we hope reach a higher level of 
precision and clarity. 

This study is the first large-scale investigation of ro- 
tation in field M-dwarfs, and provides the first useful con- 
straints on the period-mass relation for these objects after 
they have settled onto a common rotational sequence. This 
work also demonstrates the power of Kepler data and of the 
ACF method for rotation studies, and paves the way for a 
truly systematic survey of rotation rates on the main se- 
quence, from mid-F to mid-M spectral types, which we plan 
to examine in future papers. We will also seek to automate 
the ACF period verification stages for use on larger samples 
of targets, and run extensive simulations in order to quantify 
detection efficiency. 
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